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Abstract 

We propose a possible resolution for the problem of why the semicircular law is 
not observed, whilst the random matrix hypothesis describes well the fluctuation of 
energy spectra. We show in the random 2-matrix model that the interactions between 
the quantum subsystems alter the semicircular law of level density. We consider also 
other types of interactions in the chain- and star-multimatrix models. The connection 
with the Calogero-Sutherland models is briefly discussed. 
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1 Introduction 


In heavy nuclei,the complicated many-body interactions lead to statistical theories 
which explain only the average properties. One of these theories is the random matrix 
hypothesis [l||^. It supposes that the nuclear hamiltonian in a arbitrary basis of 
functions is a N x N matrix with N large and elements distributed at random. The 
joint probability function of the eigenvalues Ai,... Atv of this matrix model is given 
by: 

N 

P(Ai,... Aat) = exp(- Xj) II(Ai - Xjf (1.1) 

2=1 i<j 

where f3 = 1,2,4 for orthogonal, hermitean and,respectively unitary ensembles. In¬ 
tegrating over eigenvalues A^+i... Xn we get the joint distribution function for few 
levels: 

P(Ai,... Afc) = J dXk+i ■.. dXNPiXi, ... Xn) (1.2) 

All these joint distribution functions can be expressed in terms of the Dyson corre¬ 
lation function K{Xi, A 2 ): 

P(Ai,... Afc) = ^(-1 )"A:(Ai, A^J ... K{Xk, A^J 


where a is the permutation of k levels.In the special case k = 1 the Dyson correlation 
function coincides with level density K{X,X) = P{X). 

The density of levels for the 1-matrix model satisfies the semicircular law: 

P(A) = ^/3iV/2 - A2 

and the Dyson correlation function behaves as (cr <C A): 


K{X-]^a,X + ]^a) 


sin(7r(TP(A)) 

TTa{PN/2) 


The Dyson correlation function describes well the fluctuations of quantum sys¬ 
tems,but the semicircular law is not observed in the experimental data for the density 
of levels. A possible resolution of problem is to consider instead one random matrix 
few random matrices in interaction. As we will see , even a small interaction gives a 
calitatively new behaviour for the level density. 

As an interesting generalization of the random matrix hypothesis is to consider q 
matrices describing q nuclear systems in interaction.The total action of such system 
is: 


q-l N 

Si = Y V! Cc 

a = l 2=1 


q N 

L 

a=l 2=1 


^(«)a(“+i) 


(1.3) 


This system describes a chain of matrices with neighbour interaction.We can add a 
term describing the two-body interaction of constituent nuclear subsystems: 

N 

/ . / . Aj 

|a-/3|^li=l 
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We have different sets of energy levels ,... ,a = 1... q with distribution 
probability: 

P(aS'\ . 4^) = exp(5) n(AS'^ - A«)(aS'^) - xf) (1.4) 

i<j 

We have level repulsion only for the first and last energy level set.Hence for this 
model the intermediate energy level sets are ” classical” and interact with ” quantum” 
first and last energy level sets.Integrating over all intermediate matrices we remain 
with a two-matrix model. 

Kharchev and others have considered the so-called conformal matrix models that 
contain additional repulsion terms also for intermediate matrices ||^. 

Another special random matrix model is the star-matrix model having the action: 



2=1 a=l 2=1 a=l 2=1 


The joint distribution of this model reduces again to that of 2-matrix model. 


2 Quantum Chaos in two-matrix model 


We introduce the distribution probability: 

N 

P(Ai... Aatj/xi ... hn) = exp^(I/i(Ai) V2{ni) + CiXim) ^[(Ai - Aj)(/ri - Aij)(2T) 

2=1 i<j 


with Va{T) = -|- UaT ,« = 1, 2 and joint distribution function : 

P(Ai... Aj, /ii... fij) = j dAj+i... cZAat, ... diXNP{Xi ... Aw, /ri... ^w) (2.2) 

We show that the level densities P{X),P{fj,) and the joint probability distribu¬ 
tions P{Xi, X 2 ), P{fJ-i, fJ- 2 ) are exactly like those of the hermitean 1-matrix model with 
distribution probability (1.1): 


P{X) = Ph erm (A0,H(/r)=Pw erm 

H(Ai, A 2 ) ~ Pnerm (A^, A 2 ), PifJ'lj 1^2) ~ Ph erm (A^I , 1 ^ 2 ) 


(2.3) 


The new joint probability distributions P{X, /r) behaves in a different way because 
we have not energy repulsion between levels of different sets. 

If we set from beginning the coupling c = 0 we get two independent orthogonal 
1-matrix models and we have: 


-P(A,Ai) = Porth{^')POrthifJ-') 


For c 7 ^ 0, P{X,fi) behaves like the 1-matrix Dyson correlation function: 


P{X,^^)r.K{X,f^) 


When c —> 0,P(A,/u) does not split in two orthogonal 1-matrix models. 
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Here X', fj! are related with the coefficients of the Q-matrices. 


with; 


, _ A — ao ^ — bo 


ClU2 — 2t2Ul ClUl — 2tlU2 

“0 - -3“’ °o - 


ai = 


4tit2 - q 
2t2 


Atit 2 - cf 


M = - 


Atit 2 - cf 
2ti 


4tit2 - cf 


(2.4) 


(2.5) 


In the rest of the section we demonstrate the above relations. We introduce the 
Q-matrices,which in the two-matrix case have the form: 


Qi — 4+ + oq/o + oi€-, Q 2 — bolo + bie—, 


( 2 . 6 ) 


with 


N N 

1+ — ^ ) I^n,n+lj Iq — ^ ) En^ni 6 — 
n=0 n=0 


The Q matrices are defined as: 


N 

'y ) 

n=0 


a,nmhmj — 1? 2 


where hn = /loi?” and R = c/(4t2ti — c^)- 
T] are orthogonal polynomials 


Cn(Al) — A” 4 - ... , rjn{fJ^) — + ... 

satisfying the orthogonality condition: 

J dAdMen(A)e^l+^2+'^^^77™(/i) = hnSnm 


(2.7) 


From the definition of Q-^atrices {Qi^mn^m = Xn^n, Q 2 mn'nm = Mn?7n we have the 
following recursion relations of the orthogonal polynomials: 

^Cn(A) — (^}^^l(A) T Q'0?n(A) + Oi(^}^_l(A) 

Wn(Ai) = Vn+lilJ^) + borjnifJ^) + 6l??n-l(/U) (2.8) 

Solving these recursion relations it follows that rj are Hermite functions: 

?n(A) = CXnHn{X ),7/m(M) — ) 

To get the proportionality coefficients an, Pm we use the orthogonality relation and 
the Gauss transform: 

(27ru)-i/2 J ^ (1 _ 2ur/^Hn{{l - 2u)-^/^x) (2.9) 
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Writing the action as: 


S 


with; 


hi (A) + V2(h) + cA/U = tiA^ + ttiA + + U2fJ, + cXfj, = 

U2 + cX\‘^ / A - ao A ^ 


'S'o + ^2 ( h + 


2^2 


V V^i J 


So = 


tiU2 + t2Ui — CU 1 U 2 
Atit2 - 


( 2 . 10 ) 


( 2 . 11 ) 


we have: 

y (iA(i/rCn(A)e^i+^2+cAM^^('^^ ^ a^j3^5nme^° V4tit2 - ^ 

= ho5nmR^ = ho 4 ^ ^ 

In conclusion; 

Cn(A) = ( 27 rn!)-i/ 22 -n/ 2 (^)n^^(y)^ (2.13) 

r?^(/i) = (27rn!)-i/22—/2 (v/^)-F^(^') 


and; 

ho = (4tit2 - c^)“^/^ exp(5o) 

We can now calculate the joint probability distribution P(A, ^).Because we can 
write the two Vandermonde determinants in terms of orthogonal polynomials Vm 


A(A)A(f,) = ^ Cn(Al)Hji(A2 . . . Atv) X! ^m(Ml)0m(l^2 • • • I^^n) 

n m 

and the algebraic complements satisfy: 

r N 

/ IT (hAj(ip-j)‘^n(A2 • • • Ajv)0m(/^2 • • • /^7v) — (.A^ l)’<^nm 

i=2 

we get for joint probability distribution: 

, N-l 

P{X,^i) =—e^ ^h-^^n{X)r]n{^J) (2.14) 

n=0 

It is easy to derive the expression for symmetric joint distribution of pairs of 
eigenvalues in terms of P{X,fj,): 

i^(Ai,...,Afc,/ii,..., /ifc) = ^^(—1) P[Xi, //cti ) ■ ■ ■ P{Xk, ) 

cr 

Integrating in Aj+i ... A^ we obtain the asymmetric joint distribution of eigenvalues: 
P(Ai,..., Xj,^ii, ...,^ik)= ^(-l)'^P(Ai,/r^J ... P{Xj, 
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In the limit of large N we have the usual behaviour of semi-circular law: 


P(A) = V2N-X'^,P{n) = 

To calculate the joint distribution of two eigenvalues P{X, n) in the large N limit 
we associate it with the quantum mechanical system : 


[^(Px+pfi) + + cX^J](|)n{X)'^jJm{^J-) = Enm4>n{^)'lp7n{p) 

where px = id/dX,p^ = id/dp are the usual momenta operators and 

(f>nW = exp(-A'^/2)ry„(A) 

V’m(M) = exp(-/i'^/2)^m(Ai) (2-15) 

For c = 0 we get two decoupled quantum systems: 

{pi + a'2)0,(a) = 2f;i,„,^„(a) 

{pfi + p‘^)l/m.{p) = 2£'2,mV’m(M) (2-16) 


where Enm = -Fgn + -E' 2 ,m- 

In the large N limit Enm behaves like ~ N and because we are searching for sym¬ 
metric solutions we have Ei^n = -F 2 ,m ~ A^/2.The joint distribution of two eigenvalues 
P(A, p) will be: 


P{X,p) = ^2Ei^n - A'2^^ 


p 


/2 


(2.17) 


or 


p{x,p) = Vn - - p>^ 


(2.18) 


We can see that for c = 0, P{X, p) is the product of density energy levels for orthog¬ 
onal anssembles.If we integrate the last matrix ,we get the 1-matrix model .In our 
case this is equivalent with the condition 2172,m = + 1 ^ 2 (/^) = 0 in ( ^.16 ) or in 

other words the second system has no contribution in the joint distribution of two 
eigenvalues. The equation (2.17) is replaced by: 


P(A) = x/2N - A'2 


For c 7 ^ 0, after summing relation ( |2.14| ) and using the asymptotic formula (n 
large) for the Hermite polynomial (near origin): 

Hn = ^^^^2+^1) cos(\/2n -h 1 - ^) + 0{l/x/n) 

we obtain (up the exponent S + (A'^ -|- p'‘^)/2): 


smVm{P - p') 

P{\P) - 77777 - 7 ^’ 0 

7rN{X' — p) 


(2.19) 
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We also get for arbitrary A, //, A <C /r: 


where: 


P{X,H) ~ 


sin y/2N — (Q;A)^e(aA) 
TrNe{aX) 


e 

a 


1 1 


\f2a\ \/2h\ 

_^ 

2 \ \/2ai \/26i 


( 2 . 20 ) 


( 2 . 21 ) 


For the asymmetric potential ti = l/(a + r'f '= l/(a — r)^,(r <C a) and a 
small interaction c « 0, we have e ~ rja^^a ~ l/(2a) and e <C a. When r —> 0 
(symmetric potential) P(A, /U ~ A) tends to the level density of hermitean 1-matrix 
model PHerm{X). The interaction (even a small one) of asymmetric energy levels 
changes dramatically the level density P(A, A) of the system. 

If for T —> 0 we get the usual semicircular law, a small asymmetry creates some 
peaks in the level density P(A, A) (see figure 1). The observed behaviour is the 
quantum analog for chaotical behaviour of two interacting classical oscilators. 


3 g-matrix model 


As a random g-multimatrix model we choose the one with partition function: 


Z = / n A„A(Ai)A(Ag) exp(^ ^ c«A aAa+i) 

We show that the joint probability is : 

P{Xa, Xp) = PHerm{Xa, X'p), 1 < a < /3 < q 

where: 


(3.1) 


(3.2) 


A' 



The parameters Oq are the coefficients of the Q-matrices. 

The Q{a) have only three non-vanishing diagonal lines, the main diagonal and 
the two adjacent lines. 


Q{a) = bal+ + OaC- 


(3.3) 


where in the particular cases we know that bi = 1 and Og = R. We can write the 
parameters in terms of the determinants of two matrices (we use the results of the 
paper Q): 


R = 

d-Q — 


(-l)“(ciC 2 ...c„_i) detA„_i 
(-l)'?ciC2 .. .Cg_i(det Ag) 


(-l)"ciC2...Ca_i 


det Yr 


0 + 1 


det X„ 


(3.4) 
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The matrices and y^,are defined as follows 


and 



/ 2ti 

Cl 

0 

0 

0 

\ 


Cl 

2t2 

C2 

0 

0 


A, = 

0 

C2 

2^3 

0 

0 



0 

0 

0 

• • • 2to_i 

^k — 1 



V 0 

0 

0 

Cfc-l 

2ta 

/ 



/ 2ta 

Cfc 

0 

0 

0 


Cfc 

2ta+l 

^/c+1 

0 

0 


0 

Cfc+1 

2^q:+2 

0 

0 


0 

0 

0 

• • • 2tg_l 

Cg—1 


V 0 

0 

0 

• • • Cq—1 

2tq 


(3.5) 


(3.6) 


Of course Yi = Xg. 

As we made before for the 2-matrix model we introduce the orthogonal polyno¬ 
mials 


^n(Ai) = A” -|- lower powers, Vn{\) = A” -|- lower powers 

which satisfy the orthogonality relations 

J dXi . . . dXq^niXi)e 'Tlmi.Xq) = hnSnm 

where 

q q-l 

S - ^ ) taXiy ^ ) Cq, Aq Aq .-|_1 . 

a=l <y.=l 

We introduce also the basic intermediate functions: 

. a—1 


and 


where we denote 


Ct\Xa)= j l[dXfsU>^,)e^X 

'' 13=1 

[ n dXfse^'^-rimiXq). 
(3=a+l 


OL—l OL—l 

Sa = ^ tfsXp ^ C/3A/3A/3+1. 
13=1 (3=1 

q-l q-l 

S'a= X! + X! ^l3^l3^3+3- 

(3=a+l f3=a 


(3.7) 


(3.8) 


(3.9) 









Obviously we have 


Cn = ?n(Al), V^nH>^q) = Vm{^q)- 

In the general case when we have arbitrary potentials one sees immediately that 
and are not polynomials anymore.In our case of gaussian potentials these 
intermediate functions are again Hermite functions ,but with different arguments. 
However they still satisfy an orthogonality relation 

r ^ 

Y[ = Snmhn, I < ct < P < q. (3.10) 

These basic intermediate functions permit to write the intermediate Q matrices as: 

J = Qa,nmhm = Qa,nmhn, 1 < a < Q. (3.11) 

The equations satisfied by basic intermediate functions are: 

= Qa&\ l<a<q. (3.12) 

Aq??^“^ = 1 < a < g. (3.13) 

These equations together with the explicit form of Q-matrices permits to find the 
basic intermediate functions 

^(Aq.) — ba^n+li^a) T ®Q:Cn—1 (Aq.) 

= {aa/R)'nn+l{K) +baR'nn-l{K) (3.14) 

Solving these recursion relations it follows that , 77 ^"^ are Hermite functions for 
gaussian potentials: 

V 

Using intermediate basic functions we get for joint probability: 

. N N 0-1 

P{Xo^,X0)= {l[dx‘t^dxf^){f[ n dxt^^)x 

i=2 i=l 7=a+l 

X detKl“)(A;.“))] det^f )(Af (3.15) 

U J ij J 

Integrating over intermediate eigenvalues dx["'\'y = a + I,... P — 1 we obtain the 
joint probability of two-matrix model for which we already know the result.Hence 
we get the result ( |3.2| ). All derivation above is valid also for more general poten¬ 
tials , polynomial-like 17('r) = X)fc=i or not. The sufficient incredients are the 
coefficients of the Q-matrices. 


9 



4 Star-matrix model 

We study the star-matrix model with partition function; 


z=j iKA® n d>^) n[(\® - >f) n >< 

i=l 0=1 i<j 0=1 

N q q 

X exp(^(yo(Ar + £ F„(aS“)) + £ c^Af^Aj")) (4.1) 


2=1 


0=1 


We define the orthogonal polynomial basis as and (instead of one conjugate 
polynomial rjm q + I polynomials rjm'^: 


<? 

ri 

0=1 


dA'") n <iA'”)S(A'"))e''«+E«.<''-«-'''”''“’l n = KSn,n, 


0=1 


m = rua, a = 1,... q. (4.2) 


This basis is unusual but it works quite well at least for gaussian potentials; Vair) = 
tar'^ + UaT, a = 0,1,... g 

We introduce Q-matrices as; 


9 

n 

0=1 


dA<»' n <!A<“)S(A<»')A<“le'»+E:=.(VA+=.A<».A<.., ^ 

m = rUa, a = 1,... q. 


9 

n 

0=1 


(4.3) 


The coupling conditions are; 


qPo + 2toQo + 1*0 + E=«<3« = 0 
0=1 

P CK ~\~ ‘^tcx.Qa P ^<y. ^<y.Qo — 0, Cl — 1, . . . ^ 

With the following parametrization of Q-matrices: 

Qo = I-\-aolo-\-aiC- 

Qa — ^o/-^o-^-|- ^o-^O P -^O^—; tt — 1, . . . ^ 

we arrive at following equations: 

2taRa + Codl = 0 
2to^o “t“ U “h Cq/R(x — 0 
2t(yd(y^ U(y CotlQ - 0 

2io + E^=0 

-Tlo 

2toao + Uo + Y^ Cada = 0 

2toai + qn + ^ CaRa = 0 


(4.4) 


(4.5) 


(4.6) 
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Solving the coupling conditions we get : 


2q 1 CoiUa , 

oi = “2uo) 

h - ^ 1 R - 

Va - 0.2 V ^ I'aJi^a -. a 

-il L(y^J~L 


dry - 


1 


Atr 


{CcyUQ 2tQUct “1“ Ua ^ ^ Cq, 'y ^ 

ZLqi 


CqiUci , 

2U ' 


where A = Uq — Yj Ca I ta¬ 
in the same way we get the basic functions for Q-matrix model we can obtain 
them for star matrix model: 




(4.7) 


with: 


y(0) ^ AW-gp („) ^ AW-d^ 

Because these basic functions satisfy relation: 

fc(A'»>) = / e'A+'="A"»A'"',„(Al“)) 

we can integrate over Vandermonde determinants: 

det[r/f (Af)] = [ e'"“+'^“^‘°^^^“’det[r?i“^(A;“))] 
ij J ij 


(4.8) 


(4.9) 


(4.10) 


Then we have for the joint probability of two eigenvalues the simple expression: 

P(A(“),A(^))~PHerm(A'(“\V(^)), a,/3 = 0,l...g (4.11) 


with A',/i' given by equation ([4.81). 


5 Generalized Calogero-Sutherland model 

The connection with Calogero model permits the calculation of the joint distribution 
functions for random multimatrix models for other ensembles, different from the 
hermitean one. 

We obtain the Calogero model related to the 2-matrix model .The eigenvalue 
problem for Calogero model follows from the heat equation satisfied by the Itzykson- 
Zuber integral. 

We introduce the kernel: 

K{X,Y\t) =< X\e-^^\Y >= (27rt)-^'/2 f dUe^p[-^Tr{X - UYU+)] (5.1) 

J 
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which is related with the Itzykson-Zuber integral K{X,Y\t = 1) = exp(—+ 

F2))/(x,y): 


I{X,Y) = / d[/exp[Tr(X17yC/+)] = 


detij{e^^y^) 

(A(X)A(y))^/2 


The kernel (5.1) satisfies the heat equation 

Dx)k{X,Y\t) = 5{X,Y) 

where y|t) = (A(A)A(y))^/^Ar(A, y|t) and the laplacian is: 

9 a ^2 + 9(9 ) 


2^dx‘i 2^2 


Solving equation (^) gives: 

A(A,y|t) = ( 27 rt) 


_7V2/2 r 1 




(5.2) 


(5.3) 


(5.4) 


(5.5) 


from which follows the expression for the Itzykson-Zuber integral {a is the permuta¬ 
tion). 

We introduce the function: 


$(A|t) = J k{X,Y\t)^{Y)dY 


(5.6) 


that fulfills the heat equation with initial condition <I>(A|t = 0) = 4>(A). 

We can search for stationary solutions in the form 4>(A|t) = J2n ^n(A)e“^"^ 
where 4>„(A) satisfies the Calogero equation (without potential term): 





i<j 



^n(A) 


En^n{X) 


(5.7) 


The eigenvalues of matrix X are chosen such that yi < 1/2 ■■■< Vn- These eigenvalues 
yi ■■ .UN are mapped by the kernel K{X, Y\t) into ... x„(^x)- 

For t —> 0, the kernel k{X,Y\t) tends to “ Vi)- Hence if we 

consider 4'(A) as a particular solution of Calogero model with xi < X 2 ■ ■ ■ < xn 
,the function <I>(A|t = 0) is the general solution for eigenvalues Xi in arbitrary order, 
being the linear combination of functions ^(aX) : 

$(A|t = 0) = ^ T(aA), T(aA) = r?,T(A) 


where a is the permutation of eigenvalues x* ; r/o- = — 1 for free fermions {(3 = 2 for 
hermitean matrices) and r/o- = +1 for free bosons (/3 —> 0 for harmonic oscillator) . 

For t —> 00 the dominant contribution is given by the vacuum configuration 
<I>o(A). The kernel K{X,Y\t) plays the role of instanton propagator connecting 
the initial vacuum configuration To(y) = (A(y))^/^ to final vacuum configuration 
$o(^) = (A(A))^/2. 
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For 2-matrix model we can define the Generalized Calogero system: 


(5 (3 

-^(5^ ^ 7m2) + 9 (9 " 


2^2 


+ 


Kj 


; (Ai-Aj)2 


r)+ 


+ ^(l^l(Aj) -|- V2(^i) -|-cAi//i)J — Enm^n{^)^miM) (^’S) 

When c = 0 the generalized system splits into two Calogero systems : 

1 


--V —+ -(--l)y 
‘>^9Ar 2 4 


i<j 




i<j 


(4 - A ,)2 

1 

- i^jY 


+ ^A'2U„(A) = i?i,A(A) (5.9) 




The ground states can be written in terms of the eigenfunctions (2.14): 

$o(A) = (detijCi(Aj))^/4xp(-^ Af/2) 

i 

To(/i) = (de%r/i(^j))^/4xp(-^/rf/2) 


(5.10) 


We can see that the probability of amplitudes ( 5.10| ) is the partition function of the 
2-matrix model: 

^ ~ I “ J 


The system ( ^.10 ) permits us to calculate the joint probability P{X, p) for general 
ensemble. It coincides with formula (2.18|) (for c = 0) where we replace N by j3N/2: 


P{X,fi) = ^f3N/2 - P^^JpN/2 - ^'2 


(5.11) 


6 Conclusions 

These models present interest in the study of quantum chaos for q systems interacting 
in various ways. The density of levels depends on the total energy which behaves like 
N, for large N. The interaction of q subsystems redistribute the energy between the 
subsystems and change in non-trivial way the joint distribution functions. Different 
kinds of interaction (chain or star- type ) give different probabilities for energy levels. 
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FIGURE CAPTIONS 

Eigure 1. -Represents the level density P{x,x) in terms of the energy x = a\ and 
the asymmetry y = Ne.Fov y = 0 we have the semicircular law P{x,x) = 'f2N — x"^ 
and for small y 7 ^ 0 we get the oscilations of level density. 
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